PROCEEDINGS 

°F SCIENCE 



MKPH-T-11-21 
HIM-2011-11 
CERN-PH-TH-20 11-280 



o 

(N 

> 
O 



> 



X 



Lattice Determination of tiie Anomalous l\/lagnetic 
l\/loment of tiie l\/luon 



Michele Delia Morte, Benjamin Jager* Hartmut Wittig 

Institutfiir Kernphysik and Helmholtz Institut Mainz, 
Johannes Gutenberg-Universitdt, D-55099 Mainz, Germany 



E-mail: |norte@kph . uni-mainz . de 


jaeger@kph . uni-mainz . de 


wittigSkph . uni-mainz . de 





Andreas Juttner 

CERN, Physics Department, TH Unit, CH-1211 Geneva 23, Switzerland 



E-mail: A^ndreas . Juttner @cern . ch 



We compute the leading hadronic contribution to the anomalous magnetic moment of the muon 
flJJ^*^ using two dynamical flavours of non-perturbatively i?{a) improved Wilson fermions. By 
applying partially twisted boundary conditions we are able to improve the momentum resolution 
of the vacuum polarisation, an important ingredient for the determination of the leading hadronic 
contribution. We check systematic uncertainties by studying several ensembles, which allows us 
to discuss finite size effects and lattice artefacts. The chiral behavior of flJJ^° turns out to be 
non-trivial, especially for small pion masses. 



XXIX International Symposium on Lattice Field Theory 
July 10-16 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Lattice Determination of the Anomalous Magnetic Moment of the Muon 



Benjamin Jager 



1. Introduction 

The anomalous magnetic moment of the muon is defined as half the difference between the 
gyromagnetic factor of the muon and 2: = {g^ — 2)/2. This quantity is commonly used to 
test the Standard Model of particle physics, since it can be measured and computed to a very high 
precision. Currently a discrepancy of 3.2 standard deviations between experiment and theory is 
observed [1], which might be a hint for physics beyond the Standard Model. The hadronic contri- 
butions, especially the hadronic vacuum polarisation aj^^^, dominate the theoretical uncertainties. 
Currently the Standard Model prediction of a^^^ is determined using the optical theorem to relate 
the cross-section e^e^ — > hadrons data to the vacuum polarisation. A calculation from first prin- 
ciples is clearly desirable. Here we report on our current effort towards a precise determination of 
the leading hadronic contribution to the anomalous magnetic moment of the muon using Lattice 
QCD. 

2. Lattice Setup 





a [fm] 


lattice 


L [fm] 


m;r [MeV] 




Labels 


5.20 


0.079 


64 X 32^ 


2.5 


471 -317 


6.0-4.0 


A3 - A5 


5.30 


0.063 


64 X 32^ 


2.0 


644 - 447 


7.9-4.7 


E3 -E5 


5.30 


0.063 


96 X 48^ 


3.0 


323, 277 


5.0,4.2 


F6,F7 


5.50 


0.050 


96 X 48^ 


2.4 


541,431 


6.5,5.2 


N4,N5 



Table 1: Summary of simulations parameters. Measurements are performed on configurations separated 
by 8 units of molecular dynamics time at least. The scale and pion masses are still preliminary and taken 
from [10,11]. 

We use two dynamical flavours of non-perturbatively if?'{a) improved Wilson fermions [2] and 
include a partially quenched strange quark. Similar studies have been performed in the quenched 
approximation [3,4] and in the theory with two [5,6] and three dynamical flavours [7,8]. Our 
measurements are performed on a subset of the gauge configurations generated as part of the CLS 
project [9]. The simulations parameters are listed in Table |T| On the lattice the vacuum polarisation 
tensor is defined by a Fourier-transformation of a current-current correlator 

n^vl^) = a4£^>'?(-+''A/2-«v/2) (yc (^)_,c(o)^ ^ ^2.1) 

X 

where we use the conserved point-split current J'^ (x) [4] given in the case of Wilson fermions by 

■^mW = \(^^{^ + ai^)i^ + Yn)Ui^{x)q{x)-q{x){\-Y^)U^{x)q{x + aiX)^. (2.2) 

Preforming the Wick contractions in Equation produces connected as well as disconnected 
contributions. Disconnected diagrams are computationally expensive and neglected in this study. 
Nevertheless two-flavour chiral perturbation theory of NLO allows us to estimate the disconnected 
contributions to be —10% of the connected one [12]. 
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Current conservation implies that the vacuum polarisation tensor n^v(^) can be related to the 
vacuum polarisation Yl{q^) in the following way: 

^tiv{q) = {5^vq^ - qfiqv) n(^^). (2.3) 

For space-like momenta, the relation between the vacuum polarisation Tl{q^) and the lowest order 
hadronic contribution to the anomalous magnetic moment of the muon a^^'^ has been derived in 
[3,4,13] 

-] / dq^f{q^)n[q^), (2.4) 



■ 71 J Jo 

where the kernel of this integral is given by 

mlq^Z\l-q'Z) q^-^7^4^ 

l+mlq2Z^ ' ^ = 2^^^ ^'-^^ 

and tl{q^) = 4k (n(^^) — n(0)) . We have implemented partially twisted boundary conditions [14] 

q{x + Lk)=e'^''q{x), (2.6) 

which allow us to access any value of the momentum ^n — where H is a vector of integers. 
In this way we are able to improve the sampling with data points, in particular the kinematical 
region q^ < (^) where the kernel of the integral in Equation lA is peaked. Partially twisted 



boundary conditions can be applied to the connected diagram by reinterpreting the correlator as 
flavour non-diagonal (see [15]). 

3. Determination of a^^^ 

In order to determine the leading order hadronic contribution to via the integral in Equa- 
tion 2A, we need a continuous description of the vacuum polarisation H{q^). Therefore we perform 
correlated least square fits to our simulation data. Perturbation theory including &{as) terms [16] 
is incorporated into the fitting procedure to constrain the fit at large momentum by demanding that 
the fit-function and perturbation theory are matched at some high momentum. To ensure a smooth 
function, we also imply a matching of the first derivative at the same point. This procedure reduces 
the number of free fit coefficients by two. To evaluate the perturbative formula we use the non- 
perturbative, two-flavour Lambda-parameter parameter from [17] and the non-perturbative 
renormalization factors in [18,19,20]. We study systematic effects introduced by the choice of the 
fit procedure, by varying the fit ansatz and the matching point of perturbation theory. We chose 4 
different fit-ansatze and checked for systematic differences: 

a) a model independent Pade with degree 2 over degree 3 



b) vector dominance model including a single vector 



n{q^)=a+ (3.2) 
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c) vector dominance model with two vectors and one mass fixed to my as proposed in [7,8] 



+ 



(^2 + c2) {q2+miy 



(3.3) 



d) vector dominance model with two free vector masses 

b d 



n( 



■ a + 



+ 



(3.4) 



periodic be 
twisted be 



had 




pert 



0.2 GeV" 



Figure 1: Left: The subtracted vacuum polarisation tl{q^) computed on the F6 ensemble (/3 — 5.3; 
L = 3.0fm; m„ — 323 Me V) using twisted and periodic boundary conditions. The blue solid line shows 
the fit to the double vector ansatz c) matched to perturbation theory (light blue line). Right: The different 
contributions to from the dispersion integral broken down to several momentum ranges, shown in both 
figures indicated by different colours. 

The fits, except the vector dominance model fit, turn out to give consistent values within the 
statistical uncertainties, when the matching point of perturbation theory is chosen larger than 
1 GeV. In Figure |l| we show the subtracted vacuum polarisation together with the perturbation 
theory and fit c) (i.e. double vector). The effect of twisting is illustrated by showing periodic and 
twisted data in the same plot, which demonstrate a clear improvement on the momentum resolution 
of ti{q^). The remaining integration of Equation 2A is performed numerically. The right panel of 
Figure [I] shows the individual contributions to a^^'^ separated into different momentum regions. 
Part I shows the area for to m^, in which no data points occur. Nevertheless this small momentum 
range is constrained by the condition that ti{q^) — )• for ^ — )• and the smooth behaviour of the 
fit curve. The second momentum range (II) displays the region in which twisted data points begin 
to contribute. In the third region (III) periodic and twisted data points give a perfect description of 
the momentum behaviour. The final section (IV) shows the contribution from perturbation theory, 
which turns out to be negligible. The overall statistical error, estimated via a bootstrap procedure, 
is dominated by regions I and II and ranges from 2% to 7% for the different ensembles. 
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4. Results 

In Figure ^ we show the result for a^'^ computed for the ensembles listed in Table |l| as a 
function of m\. We observe a clear curvature and steep rise for small pion masses. To obtain 
at the physical point we need to perform a chiral extrapolation. In addition to a linear behaviour in 
m\ we also include a logarithmic term, which is motivated by chiral perturbation theory 

af^^{ml)=af+Bml + Cml\og{ml). (4.1) 

We restrict the fit to the four most chiral j3 = 5.3 data only to avoid mixing of cutoff effects and the 
chiral extrapolation. 



aJJ'-O/lO-"' JV/ = 2 aHLO/lO-i" N, = 2 + Iq 




ml [GeV^I ml [GeV^] 



Figure 2: Left: The simulation results for the hadronic contribution to a^ computed using two flavours, 
shown as function of m^. The chiral extrapolation (blue curve) is performed using an ansatz motivated by 
chiral perturbation theory. Right: Corresponding results for aj]^° including a partially quenched strange 
quark. 

We find that this extrapolation describes the whole set of data points quite well even those 
which not included in the fit. If we use this extrapolation, we obtain for a^'^ in the two-flavour 
case 

af{Nf = 2) = 537.1 (53.8)stat (37.6)chirai • 10-^°. (4.2) 
For the case of an additional quenched strange quark we end up with 

af{Nf = 2+\q)= 612.4 (49.9)s,at (48.5)chirai • 10-^°. (4.3) 

We repeat the analysis using a linear extrapolation for the 3 most chiral j3 = 5.3 data points to 
estimate the uncertainties from the chiral extrapolation by using half the difference of the central 
values. 

Since finite size effects and cutoff effects rather depend on q^, it is more instructive to study 
them using the subtracted vacuum polarisation tl{q^) than a^^. The left panel of Figure ^ shows 
the quantity tl{q^) for two different ensembles which have roughly the same pion mass and volume. 
This allows us to look for cutoff effects, which turn out to be below 3% in the momentum range 
< 1 GeV^. The right panel of Figure |3| offers an insight to finite size effects and cutoff effects, 
showing two ensemble which differ in volume and lattice spacing. It turns out that these two effects 
are around 5% and within our statistical precision. 
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Figure 3: Left: A study of cutoff effects comparing the subtracted vacuum polarisation for two ensemble 
with roughly the same pion mass m^ w 450MeV and volume L = 2.5 fm. Right: A comparison of n(^^) 
for two ensemble with the same pion mass ntji = 320 MeV. 



5. Conclusions and Outlook 

The determination of aj^^*-* using Lattice QCD is feasible, but still requires further improve- 
ments to make an impact on phenomenology. Partially twisted boundary conditions extend the 
accessible momentum range for the vacuum polarisation and thereby reduce the statistical and sys- 
tematic uncertainties. At present the individual data points for a^^'^ have statistical precision of 
2% to 7%. Summing up all individual uncertainties in quadrature we end up with an overall uncer- 
tainty of 12% for the extrapolated value at the physical point. Our study of residual systematics 
indicate that finite size effects and cutoff effects change a|^^° slightly. The value will in any case 
decrease after the inclusion of the disconnected diagrams. Further details of our study will be pub- 
lished soon [21]. In the future we will study ensembles with smaller pion masses to improve the 
extrapolation to the physical point and reduce the systematics, such as finite size effects and cutoff 
effects. Once we are confident that we can control these effects to the required level of accuracy, 
we will include a dynamical strange and charm quarks into our calculations. 
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